COMPUTING TROPICAL LINEAR SPACES 



FELIPE RINCON 

Abstract. We define and study the cyclic Bergman fan of a matroid M , which 
is a simplicial polyhedral fan supported on the tropical linear space T(M) of 
M and is amenable to computational purposes. It slightly refines the nested 
set structure on T(M), and its rays are in bijection with flats of M which 
are either cyclic flats or singletons. We give a fast algorithm for calculating 
it, making some computational applications of tropical geometry now viable. 
Our C++ implementation, called TropLi, and a tool for computing vertices of 
Newton polytopes of A-discriminants, are both available online. 



1. Introduction 

Let A be an m x n complex matrix of rank m, with columns ai, a2, . . . , a n € C m . 
We denote by M(A) its associated matroid, i.e., the matroid on the ground set 
[n] := {1,2, ... ,n} encoding the linear dependences in C m among the columns of 
A. The circuits of M(A) are then the sets C C [n] such that there is a minimal 
dependence among the columns of A of the form Yliec ^ &i ~ ®- ^ ee 0x192] for 
an introductory reference to matroid theory. 

The tropical linear space T(M) of any matroid M over the ground set [n] is the 
set of vectors v G K™ such that for any circuit C of M , the minimum min{t>i | i e C} 
is attained at least twice (i.e., there exist j,k € C distinct such that Vj — Vk — 
minjui | i 6 C}). In the case where M is the matroid associated to a complex 
matrix A € i£"ixn^ -y{M) agrees with the tropicalization of the linear subspace 
rowspace(A) C C™ (using the trivial valuation on C). We will not consider in this 
paper tropical linear spaces obtained by tropicalizing using non-trivial valuations; 
for a discussion about these more general tropical linear spaces and their beautiful 
combinatorics the reader is invited to see |Spe08| . 

Tropical linear spaces are one of the most basic objects in tropical geometry, 
and interest in them has increased substantially in the last few years. They are 
the local building blocks for abstract smooth tropical varieties, they play a key 
role in defining a well-behaved tropical intersection product, and they are central 
objects for studying realizability questions in tropical geometry (see |FR10j . [KP09J, 
[ShalOp . They are also fundamental for the study of tropicalizations of varieties 
obtained as the image of a linear subspace under a monomial map |DFS07| . It 
is thus desirable in many situations to have an explicit description of them as 
polyhedral fans, i.e., as a list of polyhedral cones in R" on which it is possible to 
perform different computations. 

There are several natural polyhedral fan structures that can be given to the 
tropical linear space of a matroid M. In [FS05 , Feichtner and Sturmfels described 
a whole family J\f of polyhedral fans, all of them supported on the tropical linear 
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space T(M). They compared these fans to the coarsest polyhedral structure on 
T(M), called the Bergman fan B(M) of M, which is induced by the normal fan 
of the matroid polytope associated to M. The finest fan structure in the family 
AT is called the fine subdivision of T(M). It was studied by Ardila and Klivans 
in |AK06| . where they used it to show that the intersection of the tropical linear 
space T(M) with the (n — l)-dimensional unit sphere is homeomorphic to the order 
complex of the lattice of flats of M, and thus to a wedge of spheres. The coarsest 
fan structure in the family Af is called the (coarsest) nested set fan of M, and 
was studied in depth in |FS05| . In particular, Feichtner and Sturmfcls proposed an 
algorithm for computing the nested set fan of M by gluing together "local" tropical 
linear spaces. In general, their algorithm has the inconveniences of having to go 
over all rank(M)! possible total orders on the elements of each basis of M, and of 
performing the computation of each maximal cone in the nested set fan a multiple 
number of times. 

In Section [2] we introduce the cyclic Bergman fan $(Af) of M, which is a 
simplicial polyhedral fan also supported on the tropical linear space T(M). The 
maximal cones of <I>(A7) are described using some interesting combinatorial objects 
that we call "compatible pairs" . We prove that the rays of $(M) are in correspon- 
dence with flats of the matroid M that are either cyclic flats or singletons, showing 
that $(M) is in general a little finer than the nested set fan of M. In Section [3] we 
present an effective algorithm for computing the cyclic Bergman fan of any matroid 
M that overcomes the difficulties present in [FS05 . We carry out a C++ imple- 
mentation of our algorithm in the case M is the matroid associated to an integer 
matrix A. The resulting software, called TropLi, computes tropical linear spaces 
with great speed. It can also be used to compute basic matroidal information about 
the matrix A, like its collection of bases, circuits, or its Tutte polynomial. TropLi 
can be obtained at the website 

http : //math . berkeley . edu/~f elipe/tropli/ . 

In Section [4] we give examples of a few computations done with it and report 
on its performance. Finally, in Section [5] we describe how our computation of 
tropical linear spaces can be used to compute vertices of Newton polytopes of A- 
discriminants. A CH — h implementation of this procedure is also available online. 

2. The Cyclic Bergman Fan 

In this section we introduce the cyclic Bergman fan <&(M) of a matroid M . It 
is a simplicial polyhedral fan supported on the tropical linear space T(M) of M 
amenable to computational purposes. 

Let M be any rank m matroid on the ground set [n] having no loops and no 
coloops. Suppose 7 C [n] is an independent set of the matroid M and e S [n] is an 
element not in 7 such that 7 U {e} is dependent. There is a unique circuit of M 
contained in 7 U {e} (containing the element e), which is called the fundamental 
circuit C(e,7) of e over 7. It can be described as 

(1) C(e, 7) = {e} U {i G 7 | 7 - {i} U {e} is independent}. 

Now, let B C [n] be a basis of the matroid M. Let C W l be the polyhedral 
cone consisting of all vectors v that make B a basis of maximal w-weight, i.e., such 
that J2ieB v i ^ s maximal among all bases of M. The set T(M)b ■— T{M) n is 
called the local tropical linear space of M around the basis B. 
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Example 2.1. Consider the 3x6 matrix 

/l 1 1 \ 
A := 1 1 1 . 
\0 -1 -1/ 

The matroid M := M(A) is in this case a graphical matroid, namely, the cycle 
matroid of the graph G presented in Figure [T] The circuits of M correspond to 
minimal cycles of G and the bases of M correspond to spanning trees of G. 



1 3 




Figure 1. A graph G 

The tropical linear space T(M) is then the set of vectors uel 6 such that V\ = v 2 , 
v 3 = v 4y and min(ui,«3, W5, Ug) is attained twice. It is naturally a polyhedral fan 
with six maximal cones, corresponding to the six posibilitics for the two positions 
where min(wi, V3, 1)5, Vq) is attained. For the basis B := {1, 5, 6}, the corresponding 
local tropical linear space T(M)b is the set of vectors v £ T(M) that satisfy W3 = 
min(i>i, V3, V5, vq), which consists of only three of the six maximal cones described 
above. Note that each of these maximal cones is in several local tropical linear 
spaces. For example, the cone described by W5 > vi — v 2 — V3 — V4 < vg is in 
the local tropical linear space corresponding to the bases {1, 5, 6}, {2, 5, 6}, {3, 5, 6}, 
and {4,5,6}. 

The following theorem and its corollary appear in the work of Feichtner and 
Sturmfels |FS05| . They show that, although the tropical linear space T(M) might 
have a complicated combinatorial structure, all local tropical linear spaces are much 
faster to compute. In order to make our paper self contained, we give here a 
completely combinatorial proof of their result. 

Theorem 2.2 ( FS05J). Let B be a basis of the matroid M. For any v £ Es, v is 
in the local tropical linear space T(M) b if and only if the minimum min{wi | i £ C} 
is attained at least twice for any fundamental circuit C over the basis B. 

Proof. Assume by contradiction that v £ £5 is such that the minimum mui{vi \ 
i £ C} is attained at least twice for all fundamental circuits C over B, but v is 
not in T(M). Let D be a circuit of M such that min{wi | i £ D} is attained only 
once, and take D containing as few elements outside of B as possible. Since D 
is not a fundamental circuit over B, the circuit D contains at least two elements 
not in B. Let a £ D be the element such that v a — minjuj | i £ D}, and let 
b £ D — B be different from a. Consider the fundamental circuit C := C(b,B) of 
b over B. Since B is a basis of maximal v- weight, Equation ([T]) implies that the 
minimum mm{vi \ i £ C} is attained at 6, that is, vi, < v c for any c £ C . In 
particular, we have that a ^ C. Applying the strong circuit elimination axiom (see 
0x192, Proposition 1.4.11]) to the circuits C and D, with the elements b £ C n D 
and a £ D — C, we get that there is a circuit E C C U D — {b} containing the 
clement a. But then E is a circuit such that min{wj | i £ E} is attained only once 
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(at i = a), and E has fewer elements outside of B than the circuit D, which is a 
contradiction. □ 



Corollary 2.3. Let B 

function f B : R m -> W 
f B (x) G K" defined by 



= {bi, • ■ • j b m } C [n] be a basis of the matroid M . The 
sending a vector x = (xx, x%, . . . , x m ) G R m to the vector 



Ub{x)\, 

is a piecewise linear homeomorphism between 
T{M) B . 



if i = bj for some j , 
ifi G [n] - B; 

J m and the local tropical linear space 



We now define the combinatorial objects that we will use to study the cyclic 
Bergman fan. Fix a basis B C [n] of M. For any k G [n] — B, denote F% := 
C(k, B) — {k} (note that Fk ^ since M has no loops). Let v be any vector in the 
local tropical linear space T~(M) b , and suppose J' is a total order on B such that 
for any a,b £ B we have v a < Vb a <j b (note that for generic v this condition 
determines J uniquely). This total order J induces a function p : [n] — B — > B 
defined by p(k) :— "J'-smallest element in Fk' . We say that p is the preference 
function induced by the total order J. According to Corollary |2.3| this preference 
function p is encoding which coordinates attain the minima described in Theorem 



2.2 



that is, min{t>i | i G C(fc, B)} — Vk — Vp(k) f° r au k G [n] — B. Let C denote the 
restriction of the total order J to the image Im(p) of p. We call the pair (p, C) a 
compatible pair (with respect to the basis B) induced by the vector v. Note that 
a non-generic vector v G T(M) b might induce several different compatible pairs 
with respect to B, corresponding to different choices of the total order J . 



Example 2.4. 



Let A be the 4x7 matrix 
/I 

A := 








2 
1 



Consider the basis B := {1,2,3,4} of the matroid M := M(A). The fundamen- 
tal circuits over B are C(5,B) = {2,4,5}, (7(6, B) = {1,2,4,6}, and (7(7, B) = 
{1,2,3,7}. Let v = (0,5,2,3,3,0,0) G R 7 . It is not hard to see that the basis B 
is a basis of maximal w-weight, so v G S_b. Since the minimum min{vi \ i G C} is 
attained at least twice for each fundamental circuit C over the basis B, Theorem 
2.2 implies that v G T(M)b- There is a unique total order J on the elements of 
B satisfying v a < vi, =>■ a <j b, namely 1 <j 3 <j 4 <j 2. The preference 
function p induced by J is then given by p(5) — 4, p(6) = 1, and p{7) = 1. The 
compatible pair (p, C) induced by v in this way consists of the preference function 
p together with the linear order 1 <£ 4. Note that a different order C on the image 
of p would not be compatible with the preference function p. In fact, if the order 
J satisfied 4 <j 1 then it would not be possible that p(6) = 1. 

Proposition 2.5. Let (p,£) be a compatible pair (with respect to the basis B). The 
set of vectors v in the local tropical linear space T(M)b that induce the pair (p, C) 
is an m-dimensional polyhedral cone T(p 1 C) C M™. Its lineality space is generated 
by the vector (1, !,...,!) G R". After modding out by this lineality space, the cone 
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T(p,C) is a simplicial polyhedral cone whose extremal rays can all be taken to be 
0/1 vectors. 

Proof. Let Q = Q(p) be the partition of the set [n] with m blocks Qb '= {b} U 
p -1 ({&}), for b £ B. Note that if a vector v £ T{M)b induces the pair (p,C) then 
v has to be constant on each of the blocks of Q, that is, for any i, j in the same block 
of Q we must have Vi — Vj. We will construct a directed caterpillar tree T = T(p, C) 
with set of vertices Q encoding all further restrictions on the coordinates of such a 
vector v: If there is a directed path in T from Qb to Qb' then v must satisfy Vi < Vj 
for i £ Qb and j £ Qv ■ 

Since £ is a total order on the elements in the image of p, it naturally induces 
a total order on the non-singleton blocks of Q. We start the construction of T 
as a directed path whose vertices are all the non-singleton blocks of Q, with their 
position in the path matching the order prescribed by C (i.e., it is possible to walk 
from Qb to Qb' if b <c V). Now, for every c £ B — Im(p), add a directed edge 
from the non-singleton block Qb to the block Q c — {c}, where b is the £-largest 
element in the image of p for which there is a k £ [n] — B such that k £ Qb (i.e. 
p(k) = b) and c £ Fk- Note that such a 6 is guaranteed to exist since the matroid 
M has no coloops. Corollary |2.3| ensures that the directed tree T constructed in 
this way encodes precisely all the conditions on the coordinates of a vector v for 
it to be a vector in the local tropical linear space T(M)b inducing the pair (p,C). 
More specifically, v £ T(M)g and v induces (p,C) if and only if v is constant on 
the blocks of Q and for any directed edge Qb — > Qb' in T we have Vb < Vy ■ 

Now, it is easy to see that the set T(p, C) of vectors v satisfying the conditions 
imposed by T is a polyhedral cone with lineality space generated by the vector 
(1, 1, . . . , 1) £ R™. Moreover, after modding out by its lineality space, the cone 
r(p, C) can be described as the positive span of m — 1 linearly independent 0/1 
vectors, as follows. Think of T as a partial order on the blocks of Q, and for any 
b £ B define Wb as the sum of all coordinate vectors such that i is in the union of 
all blocks in Q greater than or equal to Qb (according to T). The cone T(p, C) is then 
equal to the positive span of the vectors {wb \ b £ B and Wb ^ (1, 1, • • ■ , !)}■ □ 



Example 2.6. Let M and B be defined as in Example |2.4| We showed that the pair 
(p, C) is a compatible pair, where p is given by p(5) = 4, p(6) — 1, p(7) = 1, and L 



is the total order 1 <£ 4. Following the proof of Proposition 2.5 the partition Q(p) 
for this preference function is {{1, 6, 7}, {2}, {3}, {4, 5}}. The directed caterpillar 
tree T associated to the pair (p, C) is depicted in Figure[2j It encodes the conditions 
for a vector v £ M 7 for it to induce the compatible pair (p, £): v induces (p, C) if and 
only if t>3 > vi = vq = v-j < — < vi- These equalities and inequalities define 
the simplicial polyhedral cone T(p, C) in M 7 . After modding out by the lineality 
space R • (1, 1, . . . , 1), the cone T(p, C) is generated by the rays e2, e245, e$. 



We will later prove in Theorem 2.11 that the extremal rays of the cones T(p,C) 
are precisely the indicator vectors of all flats of the matroid M that are either cyclic 
fiats or singletons. 

It was pointed out to the author that the construction given in the proof of 



Proposition 2.5 of the polyhedral cones T(p, C) from the directed trees T(p, C) agrees 
with a more general construction of Postnikov, Reiner, and Williams described in 
PRW08 . In their paper, a "braid" polyhedral cone cfq C R™ is associated to 
every preposet Q on the set [n]. In the case the Hasse diagram of the preposet Q 
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Figure 2. A directed caterpillar tree T 



is a (directed) tree T, their construction of the cone oq agrees exactly with our 
construction of the cone T(p, C). 

It follows from our discussion that for any basis B of M, the local tropical linear 
space T(M) b is the union over all compatible pairs (p, C) with respect to B of the 
simplicial cones T(p,C). However, since local tropical linear spaces corresponding 
to different bases might intersect nontrivially, compatible pairs with respect to 
different bases might give rise to the same cone. In order to find a canonical pair 
representing each cone, we say that a preference function p : [n] — B — > B is 
regressive if p(k) < k for all k € [n] — B. If (p, C) is a compatible pair and p is 
a regressive preference function, we say that (p, C) is a regressive compatible 
pair. 

Theorem 2.7. The tropical linear space T(M) is the union over all bases B and all 
regressive compatible pairs (p, C) with respect to B of the simplicial cones T(p,C). 
Moreover, if(jpi,C\) and (^2,^2) ore different regressive compatible pairs (possibly 
with respect to different bases B\ and Bi), then the intersection of T(p\, Ci) and 
T(p2 1 C2) is a proper common face of both cones. 

Proof. In order to show that T{M) is the union of all cones C) with p regressive, 
let v be any vector in T(M). Let B be the first basis with respect to lexicographic 
order which has maximal w-weight (i.e., such that ^2 ieB vi is maximal). The vector 
v is then in the local tropical linear space T(M)b, so it induces some compatible 
pair {p, C) with respect to B. Assume by contradiction that p is not a regressive 
preference function, so there exists a k G [n] — B such that I :— p(k) > k. Since 
I G C(k,B), Equation ([T]) implies that B' := B — {1} U {k} is also a basis of 
M. However, we have Vk — vi, and thus B' is also a basis of maximal w-weight, 
contradicting our choice of B. 

Now, suppose that (pi,£i) and (^21^2) are any two compatible pairs. It follows 
from the description in terms of directed trees given in the proof of Proposition 



2.5 that the cones T(pi, Ci) and r(p2,£2) intersect in a common face. Moreover, if 
(pi, £1) and (f>2, £2) are distinct regressive compatible pairs then the corresponding 
directed trees T(pi,£i) and T(p2,£<2) are different. It follows that this intersection 
has to be a proper face of both cones. □ 

Definition 2.8. The cyclic Bergman fan $(M) of M is the simplicial fan in K™ 
whose maximal cones are the cones T{p 1 C) with (p,£) a (regressive) compatible 
pair. The support of $(M) is the tropical linear space T{M) of the matroid M. 
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As we will see later, the cyclic Bergman fan structure $(M) on T(M) is a little 
finer than the (coarsest) nested set structure on T(M) that was described in FS05J. 
However, working with <f>(Af) seems to be better for computational purposes, since 
its maximal cones are in one-to-one correspondence with effectively computable 
regressive compatible pairs. An explicit example of how the different fan structures 



on T{M) might look like is given in Example 2.12 



We now study the rays of the cyclic Bergman fan $(M). 

Definition 2.9. A flat F C [n] of M is called a cyclic flat if it is equal to a union 
of circuits of M. Equivalently, F is a cyclic flat if and only if F is a flat of M and 
[n] — F is a flat of the dual matroid M* . 

Lemma 2.10. Let F be a cyclic flat, and suppose I C F is an independent set 
spanning F. Then F is a union of fundamental circuits over 1 . 

Proof. Denote by U the union of all fundamental circuits over / (which are con- 
tained in F), and assume by contradiction that U C F, Since F — I C U, there 
exists some i € I such that i ^ U. Let C C F be some circuit containing i such 
that \C — I\ is as small as possible. Let a be some element in C — I, and denote 
by D the fundamental circuit of a over /. Applying the strong circuit elimination 
axiom (see [0x1921 Proposition 1.4.11]) to the circuits C and D, with the elements 
a E C D and i E C — D, we get that there is a circuit C C F containing i and 
contained in CUD - {a}, contradicting our choice of C. □ 

Theorem 2.11. The rays of the cyclic Bergman fan $(M) (after modding out by 
the lineality space generated by the vector (1,1,. ..,1) € W 1 ) are precisely the rays 
generated by the vectors ep := J2i£F e i> where F C [n] is a flat of M which is either 
a cyclic flat or a singleton. 

Proof. Recall the description of the maximal cones of $(M) and their extremal 



rays in terms of directed trees given in the proof of Proposition 2.5 If (p, C) is a 
compatible pair with respect to the basis B then from this description we see that 
all the extremal rays of the cone T(p, C) have the form M>o • ep, where F — {b} for 
some b e B or F C [n] is a union of fundamental circuits over B. Moreover, since 
ep is in the tropical linear space T(M), F must be a flat of M. 

Now, suppose F — {b} is a flat of M. If B is a basis of M containing b and v is a 
generic vector in the local tropical linear space T{M)b such that Vb = max ae s v a , 
then the singleton {b} appears as one of the leaves in the directed tree corresponding 
to the compatible pair induced by v, so ej, is an extremal ray of the corresponding 
maximal cone. 

In the case F is a cyclic flat, let B be a basis of M intersecting F in as many 
elements as possible. The vector v :— ep is then in the local tropical linear space 
T{M)b- Let J be any total order on B satisfying v a < vt a <j b, and 



let (p, C) be the compatible pair induced by J. By Lemma 2.10 the directed tree 
associated to the pair (p, C) has a node Q such that the set of elements that appear 
in nodes greater than or equal to Q is precisely F. It follows that ep is an extremal 
ray of the cone T(p, C). □ 

We now discuss how the different fan structures on T(M) that have been studied 
in the literature compare to the cyclic Bergman fan $(M). Let us assume that the 
matroid M is a connected matroid. The coarsest subdivision of T(M) is called the 
Bergman fan B(M) of M, and it was studied in |FS05| . It is the fan structure 
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on T(M) inherited from the normal fan to the matroid polytope of M. The rays 
in this fan (after modding out by the lineality space) are all the vectors of the 
form e F with F a "flacet" of M \f C [n] is a "fiacet" of M if the matroid M\F 
obtained by restricting to F and the matroid M/F obtained by contracting F are 
both connected matroids). The Bergman fan is refined by the nested set fan 
(also studied in [FS05 ), whose rays are the vectors cf with F a connected flat 
(i.e., a flat F such that M\F is connected), and whose maximal cones correspond 
to maximal nested sets of connected flats of M . This fan is in turn refined by 
the cyclic Bergman fan $(M), whose rays are the vectors with F a flat which 
is either cyclic or a singleton, and whose maximal cones correspond to regressive 
compatible pairs. Finally, the cyclic Bergman fan $(M) is subdivided by the fine 
subdivision of T(M), which was studied in |AK06j . In this fine subdivision the 
rays are the vectors ep with F any flat, and the maximal cones correspond to 
maximal chains of flats. Since the last three of these fans are simplicial fans, one 
way of measuring how different these fan structures on T(M) are is to measure how 
different the following sets of flats of M are: 

{F connected flat} C {F flat, either cyclic or singleton} C {F flat}. 

A criterion for when the Bergman fan is equal to the nested set fan can be found 
in Theorem 5.3 of [FS05] , 

Example 2.12. Let M be the graphical matroid defined in Example |2.1| The 
Bergman fan B(M) is the coarsest fan structure on the tropical linear space T(M), 
and it consists of the six maximal cones discussed in Example |2.1| This coarsest 
fan structure is also equal to the nested set fan of M. After modding out by the 
lineality space generated by the vector (1, 1, . . . , 1) g R 6 , the fan B(M) has four 
rays e±2, e^4, e$, eg € IR 6 , corresponding to the four nontrivial connected flats of 
the matroid M (which are also "flacets" of M) . There is one more cyclic flat of M 
which is not connected: the flat {1, 2, 3, 4}. This implies that the cyclic Bergman fan 
$(M) strictly refines the fan B(M). In fact, the maximal cone of B(M) described 
by vi — v 2 > v 5 = v e < t> 3 = W4 gets subdivided into two smaller cones by the new 
ray ei234 of the fan $(M). The fine subdivision of T{M) is a fan with ten rays, 
corresponding to the ten nontrivial flats of M. In the fine subdivision, each of the 
six maximal cones of B(M) gets subdivided by a new ray into two cones, to get a 
total of twelve maximal cones. 

3. Computing Compatible Pairs 

Let M be a rank m matroid on the ground set [n] having no loops and no 
coloops. The cyclic Bergman fan 3>(M) described in Section [2] allows us to develop 
an algorithm for computing the tropical linear space T(M) of M in an effective 
way. As it was discussed above, the maximal cones of 4>(M) are in bijection with 
regressive compatible pairs, so the key idea for a fast calculation of <f>(M) lies in 
coming up with a good way of computing all possible regressive compatible pairs 
with respect to a given basis B. The way compatible pairs were defined made 
use of a total order J on the elements of B to construct the pair, but it would 
not be a very good idea to go over all possible such total orders if m is not very 
small. Instead, what we do is to construct recursively each compatible pair (p, C) by 
building up p and £ at the same time. Algorithm [T] describes a general procedure 
that achieves this goal. As we mentioned before, our algorithm has two important 
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Algorithm 1: Computing the cyclic Bergman fan 3>(M) 

Input: A rank m matroid M having no loops and no coloops. 
Output: A list of all maximal cones in the cyclic Bergman fan $(M). 

for each B C [n] basis of M do 

• Compute fundamental circuits: 
for each k € [n] — B do 

|_ Compute F k := C(k, B) - {k}. 

• Compute recursively all regressive compatible pairs (p, C) with respect 
to B: 

Initialize p — and C = 0, and let k be the first element in [n] — B. 
Apply the recursive procedure Pref C) described below. 
Procedure Pref (k,p, £): 
if k — end then 
| Output the constructed pair (p, £). 
else 

if Im(p) n F k ^ then 

Define p(k) :— "^-smallest element in Im(p) n Fk" . 
Let k' by the first element in [n] — B greater than k (or 
k! = end if k is the last element in [n] — B). 
Apply Pref (k',p,C). 

for each b E F k — Im(p) smc/i i/iat 6 < k do 

for each iota/ order CJ on the set Im(p) U {6} that extends the 
total order C do 

if there is no I < k in [n] — B satisfying both b € Fi and 
b <c> p(l) then 
Define p(k) = b. 

Let k' by the first element in [n] — B greater than k (or 
k' = end if k is the last element in [n] — B). 
Apply Pref (k',p,L'). 



• Output the corresponding cones: 

for each pair (p, C) output in the previous step do 

Compute the corresponding directed tree T(p,C), as described in the 
proof of Proposition 



2.5 



Output the cone T(p, £). 



features that make it very fast compared to other existing algorithms: it does not 
have to go over all m! total orders on the elements of each basis B, and moreover, 
each cone in the fan is computed exactly once, so there is no need to store them in 
memory or compare them with previously computed cones. 

The pseudocode for Algorithm [l] deserves some explanation. Its main part con- 
sists of the recursive procedure Pref(fc,p, £), which computes for k e [n] — B all 
possible ways of defining p(k) given that we have already computed p( j) for all j < k 
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(with j € [n] — B), and that we already have a total order C on Im(p) := {p(j)}j<k- 
The block starting on Line 12 deals with the case where p(k) is defined to be an 
element already in Lm(p), in which case p(k) can only be defined as the /^-smallest 
element in Im(p) (1 -Ffc. The block starting on Line 16 deals with the case where 
p(k) is defined to be a new element not in Im(p). In this case, the condition on 
Line [l8] makes sure that the definition of p(k) will not affect the compatibility of 
the pair (p, C). 



4. TropLi: A C++ Implementation 

We developed a C++ implementation of the pseudocode described in Algorithm 
[T]for the case when the matroid M is given as the matroid associated to an m x n 
integer matrix A of rank m (having no loops and no coloops). In this case, if B 
is a basis of M and k € [n] — B, we compute the set ■= C(k, B) — {k} by first 
row-reducing the matrix A in such a way that the submatrix of A consisting of the 
columns indexed by B is the identity, and then looking at the nonzero entries in 
the column indexed by k. A few minor changes were made to the pseudocode in 
Algorithm [l] in order to improve the efficiency of our implementation. For example, 
the order of the loops described by Line [16] and Line [IT] was reversed, so that the 
amount of times the condition in Line [18] has to be checked is reduced significantly. 

In order to run through all bases B of the matroid M, our code simply lists each 
subset of [n] of size m and tests directly if the corresponding columns are a basis 
of C m . It makes use of the C++ library LEDA |Gmb) for carrying out all row 
operations on the matrix A with exact integer arithmetic. A much more effective 
way of listing all bases of the matroid M would be to make use of Avis and Fukuda's 
reverse search algorithm AF96J, which will be implemented in future versions of 
our code. 

The result is a fast software tool for computing the cyclic Bergman fan 3>(A) := 
$>(M(A)) of an integer matrix A, called TropLi. This software, together with 
documentation on how to use it, are available online at the website 

http : //math.berkeley . edu/~f elipe/tropli/ 

TropLi can also be used to compute some basic information about the matroid 
A1(A), like a list of all its bases, all its circuits, or its Tutte polynomial. 

We now present a few computations done using TropLi and report on its perfor- 
mance. All of the computations were performed on a laptop computer with a 2.0 
GHz Intel Core 2 processor and 2 GB RAM. 

Example 4.1. Let A be the 4x8 matrix whose columns correspond to the affinc 
coordinates of the 8 vertices of the three-dimensional unit cube. Running TropLi 
with this matrix A as input takes just a few milliseconds, and produces lists of all 
rays and all maximal cones in the cyclic Bergman fan <5(A). The first list shows 
that there are 20 rays in the fan $(A), each of them specified as a 0/1 vector in 
]R 8 . The second list tells us that §(A) contains 80 maximal cones, where each 
maximal cone is specified by its set of extremal rays. If instead we take A to be 
the 5 x 16 matrix whose columns are the affine coordinates of the vertices of the 
four-dimensional unit cube, TropLi still takes a fraction of a second and computes 
$(A) to be a fan with 176 rays and 2720 maximal cones. 

Running TropLi with the flag "-compare" produces a comparison between the 
cyclic Bergman fan §(A) and the Bergman fan 13(A). For the three-dimensional 
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cube these two fans are the same, and thus equal to the nested set fan. For the 
four-dimensional cube the Bergman fan has 2600 maximal cones and it is thus a 
strict coarsening of the cyclic Bergman fan, even though the cyclic Bergman fan 
and the nested set fan are still equal (see Example 5.9 in |FS05j ). Our program 
outputs a list showing which maximal cones of <&(A) are part of the same maximal 
cone in B(A). 

Example 4.2. Consider the 4 x 13 matrix 

/moo o o o o o o o o \ 

A _ I 0001 1 1 1 1 1 1 1 1 1 \ 
0100-1 -2-1 -3-2-1 I • 
\0 1 -1 -1 -2 -1 -2 -3/ 

The orthogonal complement of the rowspace of A is the rowspace of the 9 x 13 
matrix 
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The matroid M{A L ) is the dual matroid M* to the matroid M := M(A). As 
we will discuss in Section [5] computing the tropical linear space of this dual ma- 
troid can be used as the key ingredient for computing the tropicalization of the 
A-discriminantal variety. For this matrix A, this variety is the hypersurface defined 
by the condition on the coefficients of a general affine linear form l(x,y) and a 
general cubic polynomial g(x,y) so that the curves l(x,y) — and g(x , y^ 1 ) = 
are tangent. 

A Maple implementation of the algorithm described in FS05J for computing 
tropical linear spaces locally takes many hours to compute T{M*). As mentioned 
above, for each basis of the matroid M* (there are 430 of them) it has to go 
through all 9! = 362 880 possible ordcrings of the rows of A- 1 . It also computes 
each maximal cone several times, so it has to compare each cone produced with the 
list of previously computed cones to see if it is a new cone or not. Running TropLi 
on the matrix A 1 - computes the 9 dimensional fan $(Af*) (which has 29 rays and 
2466 maximal cones) in less than a second. 

Example 4.3. Let A be the 5 x 20 matrix 

(11111111110000000000\ 
00000000001111111111 \ 
02001001100200100110 . 
00200101010020010101 / 
00020010110002001011/ 

The orthogonal complement of the rowspace of A can be described as the rowspace 
of the 15 x 20 matrix 
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Again, computing the tropicalization of the rowspace of ^4 can be used for studying 
the tropicalization of the A-discriminantal variety, which now corresponds to the 
condition on the coefficients of two general quadratic polynomials in three variables 
for their corresponding surfaces to be tangent. 

Running TropLi with the matrix A as input computes the 15-dimensional fan 
&(A - L ) C M 20 , which has 172 rays and 475 722 maximal cones. All the computation 
takes just 60 seconds. A lot of this time is actually spent testing which subsets of 
the set {1, 2, . . . , 20} of size 15 are bases of the matroid M(A ± ) by row-reducing 
all 15 x 15 submatrices of A- 1 . However, since M(A ± ) is the dual matroid to 
M (A) , we can instead use A to compute the bases of M (by reducing the 5x5 
submatrices of A) and then take their complement to get the bases of M{A- L ). Also, 
the fundamental circuits Cjvf*(fc, B) in M* = M(A ) can be computed from the 
fundamental circuits in M, since j £ Cm* (k, B) if and only if k £ Cm{j, [ft] — B). 
This method of computing <f>(A ) has been implemented in TropLi, and can be 
accessed by running it on the matrix A using the flag "-dual" . In this way, the 
computation of ^(A- 1 ) is even shorter: 30 seconds. 

Example 4.4. Let A be the 6 x 30 matrix 



111111111100000000000000000000\ 
000000000011111111110000000000 \ 
000000000000000000001111111111 
020010011002001001100200100110 
002001010100200101010020010101 / 
000200101100020010110002001011/ 



The orthogonal complement of the rowspace of A is then the rowspace of a 24 x 
30 matrix A . The A-discriminant is in this case the condition on the Cayley 
octad obtained as the intersection of three general quadratic surfaces in 3-space 
(see |PSV11| ) to acquire a double point. Running TropLi on the matrix A we see 
that $(A ± ) is a 24-dimensional fan in K 30 having 929 rays and 154495 683 maximal 
cones. The computation of all these 150 million cones takes a little more than 5 
hours. If instead we run TropLi with the matrix A as input and using the flag 
"-dual", the computation of $(A ) takes less than 4 hours. 



Let A be an m x n integer matrix of rank m with columns a 1: . . . , a n £ Z m , and 
suppose that the vector (1,1,..., 1) is in the rowspace of A. The columns of A de- 
termine a collection of Laurent monomials x ai , . . . , x a?l in the ring Clxf 1 ,..., x^] 
in a natural way. Consider the space <C A of all Laurent polynomials whose support is 
contained in this set of monomials, i.e. , polynomials of the form /(x) = y^ ._ 1 Ci -x ai , 
where the CjS are complex coefficients. Let C C A be the Zariski closure of the 
set of all / in C A that define a singular hypersurface in the torus (C*) m , that is, 
for which there exists z £ (C*)" 1 such that 



The variety V a is an irreducible variety defined over Q, called the A-discrimi- 
nantal variety. When V^i C C A is a subvariety of codimension 1, the irreducible 
integral polynomial in the coefficients of / that defines is called the A- 
discriminant (A^ is defined up to sign). 

Of special interest is the case when A is a general matrix whose first s rows 
ri, . . . , r s £ r L n are given by Tj = Y^iei e * ^ or some partition . . . , I s } of [n] (see 



5. An Application: Computing A-Discriminants 



df_ 

dxi 



(z) = for all i = 1, 2, ... , m. 
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examples 4.2 4.3 and 4.4). In this case, the space C consists of polynomials of 
the form 

/ = xi ■ fi(x s+1 , . . .,x m ) H hi s - f s (x s+ i,. . . ,x m ), 

where fj is a polynomial on the variables x' = {x s+ i, . . . , x m } whose support 
is contained in the set of monomials determined by the submatrix of A with rows 
indexed by {s + 1, . . . , m} and columns indexed by Ij . The A-discriminantal variety 
is then the Zariski closure of the set of such polynomials /i(x'), . . . , / s (x') that 
have a common root in the torus where their gradient vectors {dfj/dx) x ^ are 
linearly dependent. In the case where s = 2 this corresponds to the condition 
on the polynomials /i,/2 for their corresponding hypersurfaces to be tangent. If 
s = |x'| then this is the condition on the polynomials /i, ••.,/« for the (finite) 
variety that they define to have a double point. In this case, the A-discriminant is 
also called their mixed discriminant (see [CCD+p . If s > |x'| + 1 and the matrix 



A is essential (see DFS07J), then this is simply the condition on the polynomials 
/i, . . . , f s for them to have a common root, so the A-discriminant is the same as 
their resultant. An extensive geometric treatment of all these notions can be found 
in |GKZ08j . 

Computing A-discriminants is in general a very hard computational task. Even 
for very small matrices A, the degree of and its number of monomials can be 
quite large. From the definition, A-discriminants can in principle be computed by 
solving an elimination problem in the ring Cfx^ 1 , . . . , x^], but due to the huge 
size of these polynomials a Grobner bases approach does not go too far. 

In [DFS07 , Dickenstein, Feichtner, and Sturmfels proposed a way of getting a 
handle on ^4-discriminants via tropical geometry. They proved that if A 1 - denotes 
a Gale dual of the matrix A, i.e., an (n — m) x n matrix whose rowspace is equal to 
the orthogonal complement of the rowspace of A, then the tropicalization T(Va) 
of the variety Va can be computed as the Minkowski sum of the tropical linear 
space T(M(A ± )) and the rowspace of A. In the case where V a has codimension 
1, this tropicalization T(V a) is equal to the (n — 1) -dimensional skeleton of the 
normal fan of the Newton polytope NP(Aa) of A^. They used this to describe a 
"ray shooting" algorithm to recover vertices of NP(Aa) from T(Va), which goes 
as follows (see Theorem 1.2 in |DFS07j ). Assume the columns of A span the integer 
lattice Z m . Suppose w is a generic vector in R™, and let u £ Z™ be the vertex of 
NP(Aa) minimizing the dot product u ■ w. Then u can be computed as 

(2) m = ^ \det{A t ,a 1 , . . . ,<r n _ m _i,ei)| , 

where Ci_ w denotes the set of all maximal cones a in the nested set fan of M(A ± ) 
satisfying (w + M>o ■ e^) n (a + rowspace A) ^ 0, and <7i, . . . , er„- m -i are the 0/1 
extremal rays of the cone a after modding out by its lineality space. An essential 
component in this procedure for computing vertices of NP(Aa) is to compute the 
tropical linear space T(M(A' L )). It is possible, however, to replace in this formula 
the nested set fan of M(A ± ) by the cyclic Bergman fan ^(A- 1 ), which can be 
computed more easily. 

Based on our implementation TropLi for computing cyclic Bergman fans, we de- 
veloped a C++ code that computes vertices of Newton polytopes of ^4-discriminants 
in the way described above. Given an integer matrix A and a vector w € Z™, it 
computes a vertex u of NP(Aa) minimizing the dot product u-w. In the case where 
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w is not generic and this minimum is attained at several vertices of NP(Aa), the 
code uses a symbolic perturbation approach to compute one of these vertices at 
random. This software tool can also be obtained at the website 

http : / /math . berkeley . edu/~f elipe/tropli/ . 
Example 5.1. Consider the 4 x 16 matrix 

/11111100 \ 

4 _ 00 1 1 1 1 1 1 1 1 1 1 \ 

A — I 1 1 2 -1 -2-1 -3-2-1 I 

\0 1 2 1 -1 -1 -2 -1 -2 -3/ 

Running our program with the matrix A as input and the flag "-random 100" 
computes 100 random vertices of the Newton polytope of the A-discriminant A^. 
In this case, Aa is the condition on a general quadratic polynomial f(x,y) and a 
general cubic polynomial g(x,y) for the curves f(x,y) — and g{x^ 1 = to 

be tangent. Our code computes the fan ^(A- 1 ) in the same way TropLi does (with 
the matrix A as input and using the flag "-dual"), and for each maximal cone a 
computed it checks if a + rowspace(A) C R 16 has codimension 1. If this is not 
the case then the cone a will not contribute to the sum in Equation [2] The cyclic 
Bergman fan has 18 045 maximal cones, 6 675 of which have codimension 1 

after adding the rowspace of A. This initial computation takes 22 seconds. The code 
then performs the ray shooting algorithm for 100 random values of the vector w 
using the 6 675 cones computed before, and outputs the corresponding 100 vertices 
u of the Newton polytope of A^. It also prints the A-degree of the A-discriminant 
Aa, i.e., the vector A ■ u for u any point of NP(Aa) (this does not depend on the 
choice of u). In this case, the A-degree is equal to (24, 22, —6, —6), so in particular 
we see that A^ is a homogeneous polynomial of degree 46. This second part of the 
computation takes 7 minutes. 

Of course one would like to compute all the vertices of the Newton polytope 
NP(Aa)- In general, however, due to the very large number of vertices of NP(Aa) 
and number of maximal cones in T(Va) ? one has to be quite clever about the way 
the vectors w are chosen. Choosing them at random and waiting until all vertices of 
NP(Aa) have been computed is in general not viable. A great example illustrating 
all these difficulties and a few ways to overcome them is given in |CTY10) . 

Very recently an effective algorithm for recovering the normal fan of the Newton 
polytope of a polynomial from the support of its tropical hypersurface was proposed 
in |JY11) . This algorithm has already been implemented in the software package 
Gfan | Jen) . It can be used to take the description of T(Va) as a sum of a tropical 
linear space and a classical linear space and compute from it the normal fan of the 
Newton polytope of the A-discriminant A^. From this normal fan it is possible 
to recover the exact coordinates of the vertices of the Newton polytope of Aa by 
keeping track of the multiplicities of the codimension 1 cones. 

Now, even if we have a list of all the vertices of the Newton polytope NP(Aa), 
recovering the polynomial A^ is no easy task. One way to do it is to consider a 
generic polynomial whose monomials correspond to the lattice points in NP(Aa), 
and then imposing the condition that it vanishes on the image of the rational 
parametrization of Va given in Proposition 4.1 of DFS07]. This translates into 
a linear system of equations on the coefficients of this generic polynomial whose 
solution space corresponds to the coefficients of the A-discriminant A^. Note that, 
however, this procedure requires first computing all lattice points in the polytope 
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NP(Aa) and then solving a very large system of linear equations, so a more effective 
approach would be desirable. 
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